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We create an ultracold dense quantum gas of ground state molecules 
bound by more than 1000 wavenumbers by stimulated two-photon 
transfer of molecules eissociated on a Feshbach resonance from a 
Bose-Einstein condensate of cesium atoms. The transfer efficiency 
exceeds 80%. In the process, the initial loose, long-range electro- 
static bond of the Feshbach molecule is coherently transformed into 
a tight chemical bond. We demonstrate coherence of the transfer in 
a Ramsey-type experiment and show that the molecular sample is 



not heated during the transfer. Our results show that the prepara- 
tion of a quantum gas of molecules in arbitrary rovibrational states 
is possible and that the creation of a Bose-Einstein condensate of 
molecules in their rovibronic ground state is within reach. 

Ultracold samples of molecules are ideally suited for fundamental studies in physics 
and chemistry, ranging from few-body collisional physics [H El [3l E] , ultracold chemistry 
[5], high resolution spectroscopy ^ (7j, to quantum gases and molecular Bose-Einstein 
condensation [8] and quantum processing pjj. For many of the proposed experiments full 
control over the molecular wave function in specific deeply bound rovibrational states 
is needed. High densities are required for molecular quantum gas studies. Only in the 
rovibronic ground state, i.e. the lowest energy level of the electronic ground state, is 
collisional stability assured. However, direct molecular cooling towards high phase space 
densities seems yet out of reach pTQj, whereas techniques like Feshbach association [Tl] 
and photoassociation [12j either only produce molecules in weakly bound rovibrational 
levels or suffer from low production rates and low state selectivity. In order to produce 
a quantum gas of molecules in their absolute ground state, Jaksch et al. [13j proposed 
a scheme for homonuclear alkali molecules in which the technique of stimulated two- 
photon transfer is repeatedly applied to molecules associated from a high-density sample 
of ultracold atoms. The initially very loosely bound molecules are transferred in successive 
steps to the rovibrational ground state of the singlet X^'^^ molecular potential. The 
advantage of this scheme is that it is fully coherent, not relying on spontaneous processes, 
and that it involves only a very small number of intermediate levels. It promises that 
a ground state binding energy of typically 0.5 eV can be carried away without heating 
the molecular sample. It essentially preserves phase space density, allowing the molecular 
sample to inherit the high initial phase space density from the atomic sample. However, 

2 



to realize this scheme, several challenges have to be met. First, there is a large difference 
in internuclear separation that has to be bridged, i.e. the overlap between the radial wave 
function of the least bound molecules with the radial wave functions of deeply bound 
molecular levels is extremely low, potentially leading to prohibitively low transition rates 
for the two-photon transitions. Second, the scheme requires the identification of suitable 
intermediate molecular levels while strictly avoiding parasitic excitations. Third, a large 
difference in binding energy has to be overcome. On a more technical side, the lasers 
driving the two-photon transitions at widely different wavelengths need to have extremely 
low relative short term phase jitter and high long term frequency stability to allow for 
coherence and reproducibility. In important experiments, Winkler et al. [H] and recently 
Ospelkaus et al. [T3] demonstrated highly efficient two-photon transfer into lower lying 
molecular levels starting from weakly bound dimer molecules, which were associated from 
ultracold atoms on a Feshbach resonance [TT]. However, the transferred molecules are 
still weakly bound. Their binding energy is less than 10~^ of the binding energy of the 
rovibrational ground state, and wave function overlap with this state is still negligible. 

Here we demonstrate the crucial step towards full control of the molecular wave func- 
tion and towards the formation of a Bose-Einstein condensate (BEG) of molecules in their 
rovibronic ground state by linking weakly bound molecular states with deeply bound rovi- 
brational states. We coherently transfer an ultracold quantum gas of weakly bound cesium 
Feshbach molecules to the rovibrational level |z/ = 73, J = 2> of the singlet X^'^^ potential, 
bound by 1061 cm~^, or /ix 31.81 THz, more than one fourth of the binding energy of 
the rovibrational ground state. As important ingredients, we overcome low wave function 
overlap by using a suitable intermediate excited molecular state while avoiding excitation 
into loss channels, and we reference the transfer lasers to a frequency comb, allowing us to 
flexibly bridge binding energy differences of more than 1000 cm~^ while assuring relative 
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coherence. 

Fig. 1 shows the energy of the relevant molecular and atomic states. Our experi- 
ment starts with a cigar-shaped BEG of cesium atoms in the lowest hyperfine sublevel 
F = 3,mF = 3 in an optical dipole trap. For BEG production, we essentially follow the 
procedure detailed in Ref. [IS]. For Feshbach molecule production out of the BEG, we 
ramp up the offset magnetic field from the initial value of 2.1 mT to about 5.0 mT in 
10 ms. We then ramp down, sweeping across a d-wave Feshbach resonance at 4.8 mT 
after about 1 ms as shown in Fig. IB [T71 [TH]. Our procedure [IT] gives an ultracold 
and dense sample of up to 11000 molecules every 10 s. Upon lowering the magnetic field, 
the molecules are transferred from the initial state \d> to a still weakly bound s-wave 
molecular state \s> of the lowest hyperfine channel {Fi = 3, F2 = 3) via an avoided crossing 
|18] . The index i = l,2 denotes the i-th atom. Upon further lowering the magnetic field 
to about 2.2 mT, the molecules enter into a closed channel s-wave molecular state \a> 
via a second, broad avoided crossing [18]. This state belongs to the uppermost hyperfine 
channel {Fi = 4, F2 = 4) and thus has an effective binding energy of more than 2 x hucs- 
Here h is Planck's constant and z/cs ~ 9.19 GHz is the Gs clock frequency. Similar to 
\s> this state is a mixture of the ground state and the lowest triplet a^S+ state, 

coupled by hyperfine interaction, and it has zero rotational angular momentum. At a 
field of 1.9 mT, it has a binding energy of 5 MHzx/i with respect to the -F = 3, mp = 3 
two-atom asymptote [H]. As one might expect, we find that optical transition rates as 
measured below are improved when using this effectively more deeply bound state as the 
initial state for two-photon transfer instead of state |s>. We shut off the trap and perform 
all subsequent experiments in free flight to reduce the particle density in order to avoid 
collisional processes. We detect molecules in \a> via states \s> and \d> by reversing 
the magnetic fleld ramp, dissociating them on the Feshbach resonance at 4.8 mT, and 



4 



imaging the resulting atoms [T7] . 

Efficient two-photon transfer via the stimulated Raman adiabatic passage (STIRAP) 
technique [201 [H] relies on a suitable choice for the excited state |e>. In our case this state 
must have singlet character so that it can be used as a transfer state to deeply bound 
levels of the -^^S^ potential. In general, it must be well separated from other states, 
which otherwise could be off-resonant ly excited. It should thus be situated far to the red 
of the excited Si+Pi potential asymptote to avoid the high density of excited molecular 
states near that asymptote. We have performed optical loss spectroscopy starting from 
state \a> in the wavelength range from 1120 to 1130 nm, about 2300 cm~^ to the red of 
the cesium Di line. For this we recorded the number of remaining molecules in \a> as 
a function of excitation wavelength and found two progressions of lines, which we assign 
to the potential curves of the mixed — b^Hu) 0^ excited states and to the (1)^S+ 

excited state, respectively. For the present experiments, we choose for |e> a level of the 
Oj progression which is 8879.63(1) cm~^ above the F = 3,mi7 = 3 two-atom asymptote, 
corresponding to a transition wavelength of 1126.173(1) nm, see Fig. lA. We measure all 
wavelengths on a home-built wavemeter. We identify this previously unknown level as 
the 225th one of the Oj system, with about two units of uncertainty. 

The ground state level \g> with vibrational quantum number z/ = 73 is well known 
from conventional molecular spectroscopy [211 [22] ■ However, its binding energy, as well 
as the binding energy of all vibrational levels, has only been known with an uncertainty 
of about ±0.45 cm"^ prior to the present experiments [22] • We search for \g> by exciting 
the transitions from \a> to |e> with laser Li and from |e> to \g> with laser L2 simul- 
taneously [23]. The two light fields create a molecule- molecule dark state. The molecules 
initially in \a> are lost unless the second laser L2 is on two-photon resonance, provided 
that the Rabi frequency Q2 on the second transition is equal to or greater than Qi, the 
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Rabi frequency on the first transition. For coherence, stability, and reproducibihty, we 
lock both lasers to independent narrow-band optical resonators, which we reference to an 
optical frequency comb. The comb is not calibrated, but it allows precise differential fre- 
quency measurements and provides long-term stability needed for systematic line searches 
jl9j . We find the resonance condition with vibrational level z/ = 73 at 1005.976(1) and 
1005.982(1) nm, corresponding to rotational quantum numbers J = and 2. Identification 
of J is possible since the rotational energy splitting is well known. Fig. 2A and B show 
typical molecular dark resonances when we set L2 on resonance and step the detuning 
Ai of Li near 1126.173 nm. Fig. 2C shows a dark resonance involving v = 73, J = 2 
using a different excited molecular state |e'>, which is excited with Li near 1123.104 nm. 
Fig. 2D-F show dark resonances involving the neighboring vibrational levels i> = 71 and 
u = 72. These X^S+-levels were easily found on the grounds of previous data [22]. We 
determine the binding energy of these levels with respect to the atomic Fi = 3, -F2 = 3 
asymptote at zero magnetic field to be 1060.9694(10), 1088.3101(10), 1115.9148(10) cm^^ 
for z/ = 73, 72, 71 with J = 0, respectively. The binding energy of the rovibrational ground 
state z/ = is thus 3628.7053(10) cm~^, which represents an improvement of 3 orders of 
magnitude compared to the previous determination [22]. Fitting the data for the dark 
resonances with a three-level model taking into account off-resonant excitations and laser 
line widths, we determine the molecular transition strengths as given by the normal- 
ized Rabi frequencies for the transitions \a> to |e> and |e> to |z/ = 73, J = 2> to be 
fii = 27r X 2 kHz y^//(mW/cm^) and ^2 = 27r x 11 kHz ^//(mW/cm^), respectively. A 
comparison with a typical atomic transition strength of = 27r x 5 MHz ^//(mW/cm^) 
giving < 10^® reflects the minuteness of the wave function overlap. 

We are now in a position to carry out coherent transfer using the STIRAP technique. 
For \g> we choose the vibrational level with z/ = 73, J = 2. This level will allow us to reach 
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the rovibrational ground state u = 0, J = with a second STIRAP step in view of the 
selection rule AJ = 0, ±2. STIRAP uses a counterintuitive overlapping pulse sequence in 
which first L2 is pulsed on and then Li. As is well known [20j, STIRAP relies on the 
existence of a dark state of the form \D>= a{t)\a> +j3{t)\g>. With sufficient adiabaticity, 
the function decreases smoothly from 1 to 0, while the function increases 

smoothly from to 1. The initial state |a> is thus rotated via \D> into the final state 
\g>. The criterion for adiabaticity reads TpQ"^ ^ (27r)^r, where Tp is the pulse overlap 
time, ^ ^ is the peak Rabi frequency during the pulse, and F ^ 27r x 4 MHz 
is the (spontaneous) decay rate from the upper state |e> as determined from our loss 
measurements. This criterion is quite stringent, in particular in view of the low wave 
function overlap that enters into Q, while an upper (experimental) limit for Tp is given 
by the relative laser coherence time for Li and L2. We choose Tp to be approximately 
10 /is. To get sufficiently high Rabi frequencies given a maximum laser power of up to 26 
mW for Li and 5 mW for L2 we focus both laser beams to a 1/ e^-waist of about 25 fim. 
For detection, we apply the reverse STIRAP sequence after a waiting time t^ ~ 10 fis to 
transfer the molecules back into \a>. During this time we leave laser Li on to assure that 
all possible residual population in state \a> is removed. 

We perform double STIRAP about 3 ms after the production of the Feshbach molecules. 
Fig. 3A shows the molecular population in |a> as a function of the STIRAP time r, and 
Fig. 3B shows the timing sequence for the double transfer scheme. For recording the time 
evolution of the population we interrupt the transfer process after time r and measure 
the remaining population in \a>. The molecules in \a> initially disappear during the first 
STIRAP sequence. They are now in level |z/ = 73, J = 2> of the singlet X^H'^ potential. 
Then a large fraction of them returns in the course of the reverse STIRAP sequence. For 
this particular measurement both lasers are on resonance. The peak Rabi frequencies are 
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f2i 27r X 3 MHz and ^2 ~ Svr x 6 MHz. We typically obtain an overall efficiency of more 
than 65% for the double transfer process, corresponding to single pass efficiencies of more 
than 80%, assuming equal efficiencies for both passes. Fig. 3C shows the double pass 
efficiency as a function of detuning A2 of laser L2. Simulations for the three-level system 
show that the FWHM-width of the efficiency curve of about 800 kHz is compatible with 
a combination of laser power broadening and Fourier broadening. Our simulations also 
show that higher transfer efficiencies can be expected for an optimized STIRAP pulse 
sequence in which both peak Rabi frequencies are equal. 

We demonstrate coherence of the transfer process in a Ramsey-type experiment |14j . 
For this, we halt the transfer process during the first STIRAP sequence after r ^ 12 /is 
when a balanced superposition of \a> and \g> has been created with |a(T)p ~ | ~ 
by shutting off both lasers instantaneously. After a hold time Th we resume the STIRAP 
transfer, with the role of lasers Li and L2 reversed. Thus, for = the population will 
simply be rotated back into the initial state. A three-level calculation shows that the 
population in the initial state \a> is expected to oscillate at the rate of the two-photon 
detuning IA2 — Ai|/(27r). Fig. 4A shows the initial state population for Ai ^ and 
A2~27rx 113 kHz as a function of Th- Indeed, the population oscillates with a frequency 
at IA2 — Ai|/(27r), however with marked increase in phase jitter on the time scale of 30 /is. 
We attribute this apparent loss of phase coherence to a slow relative frequency drift of 
lasers Li and L2, leading to a slightly different two-photon detuning from shot to shot. 
In Fig. 4A, we have added a region indicating a frequency jitter of ±6 kHz. This value 
is compatible with the present stability of our lasers. Note that the frequency drift does 
not affect an individual STIRAP process as the transfer efficiency is very robust against 
laser detuning as shown in Fig. 3C. 

We now show that the molecular sample is not heated during the transfer process 
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and is indeed in the quantum gas regime. For this we measure and compare the rate 
of expansion of the molecular sample in state \a> without and with the double transfer 
process. In our regime the energy scale for expansion is usually set by the mean field of 
the BEC, resulting in typical expansion energies for the atoms in the range from /c^ x 2 
nK to ks X 10 nK, where ks is Boltzmann's constant, depending on the strength of the 
atomic interaction [24]. We find that the initial magnetic field ramping excites collective 
motion of the BEC in the form of a breathing mode as a result of a change in the atomic 
interaction strength (TB] . The breathing is transformed into expansion of the sample when 
the trap is shut off. We follow the expansion by monitoring the change of the Thomas- 
Fermi radius r of the sample. Fig. 4B shows this radius along the horizontal direction 
as a function of expansion time without and with STIRAP. Without STIRAP, we obtain 
from a linear fit an expansion rate of dr /dt = 1.0(1) mm/s, corresponding to an energy of 
kB X 14(4) nK. With STIRAP, the rate is dr /dt = 0.7{l) mm/s, corresponding to an energy 
of ks X 7(2) nK. Both values are compatible with the expansion of the atomic sample for 
the same magnetic field ramp. Interestingly, the rate for the case with STIRAP is lower. 
We speculate that STIRAP with the tightly focused laser beams Li and L2 preferentially 
transfers molecules in the center of the sample and is hence responsible for some selection 
in velocity space. 

It is now possible to add a second STIRAP step for transfer into the rovibrational 
ground state z/ = 0, J = 0. A suitable two-photon transition at readily available laser wave- 
lengths is via the 63rd excited state level of the 0^ potential near 1339 nm (up) and 999 
nm (down) with comparatively good wave function overlap at the level of ^ 10~^ 

and 10~^, respectively. We expect that initial spectroscopy near 1339 nm and subsequent 
search for dark resonances will be straightforward as now all transition energies are known 
to 10"'^ cm^^. Molecules in u = 0, J = cannot further decay into a lower state upon a 
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two-body collision, and they are thus expected to allow the formation of an intrinsically 
stable molecular BEC. The high speed of our STIRAP transfer will allow us to per- 
form in-situ as well as time-of-flight imaging for direct characterization of the spatial and 
momentum distribution of the molecular ensemble. With our technique any low-lying 
vibrational state can be coherently populated in a controlled fashion with full control 
over the rotational quantum number, allowing, e.g., state-specific collisional studies and 
high-precision molecular spectroscopy with possible implications for fundamental physics 
[HI E] . Our procedure can be adapted to other species, in particular to heteronuclear alkali 
dimers such as RbCs [2^ and KRb [15] for the creation of dipolar quantum gases |26j . 
For heteronuclear alkali dimers a single two-photon transfer step might suffice as a result 
of favorable wave function overlap . We expect that the combination of our technique 
with Feshbach molecule production out of a Mott-insulator state in a three-dimensional 
lattice [2B] will increase the initial Feshbach molecule production efficiency, avoiding col- 
lective excitations as a result of magnetic field ramping and inhibiting collisional loss, and 
will provide full control over all internal and external quantum degrees of freedom of the 
ground state molecules. 
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Fig. 1. A Molecular level scheme for Cs2. Molecules in a weakly bound Feshbach level 
are transferred to rovibrational level |z/ = 73, J = 2> of the singlet potential with a 

binding energy of 1061 cm~^ in a two-photon STIRAP process with wavelengths near 1126 
nm and 1006 nm via the 225th level of the electronically excited — b^Hu) 0^ poten- 

tials. The A^S+ potential has about 155 vibrational levels. B Zeeman diagram showing 
the energy of all relevant weakly bound molecular levels for initial Feshbach molecular 
state preparation [18]. The binding energy is given with respect to the -F = 3, mi? = 3 two- 
atom asymptote. The molecules are produced on a d-wave Feshbach resonance at 4.8 mT 
(see inset) and then transferred to the weakly bound s-wave state \s> on an avoided state 
crossing. Further lowering of the magnetic offset field to 1.9 mT transfers the molecules 
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from |s> to state |a>, the starting state for the STIRAP transfer. C STIRAP transfer 
scheme [20]. The molecules are transferred from the initial state \a> to the final state 
|(7>= |z/ = 73, J = 2> by means of two overlapping laser pulses for which laser L2 is pulsed 
on first and then Li. The detunings and Rabi frequencies of Li are Aj and Qi, i = 1,2. 
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Fig. 2. Dark resonances for vibrational levels v = 71,72, and 73. Laser L2 is held on 
resonance, while the detuning Ai of Li is scanned. We record the number of molecules in 
\a> while both lasers are pulsed on simultaneously. A, B, and C show dark resonances 
involving z/ = 73 for excitation with Li near 1126 nm into J = and 2 and for excitation 
with Li near 1123 nm into J = 2, respectively. D, E, and F show the neighboring levels 
u = 71 and 72 for excitation near 1123 nm. The solid line in B is the result of a three- 



level model calculation matched to the data giving Qi = 27Tx2 kHz y/i/(mW/cm^) and 
fl2 = 27r X 11 kHz \J I2/ {vsiW / cm') for a pulse time of 5 /is at intensities of Ji = 4 x 10^ 
mW/cm^ for Li and I2 = 2x 10^ mW/cm^ for L2 assuming a laser line width of 2 kHz. 
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Fig. 3. STIRAP transfer from the weakly bound state \a> to the deeply bound state 
\g>= |z/ = 73, J = 2> and back to \a>. A shows the number of molecules in state \a> as 
a function of STIRAP time r for Ai ^ ~ A2. The measured pulse overlap begins at 
5 /is and ends at about 15 /is. The second pulse overlap starts at 25 /is and ends at about 
33 fis. B schematically shows the timing for the Rabi frequencies Qi, i = 1,2, during the 
double STIRAP sequence. Laser Li is left on after the first STIRAP sequence to clear 
out any remaining population in \a>. C Double STIRAP efficiency as a function of the 
detuning A2 of laser L2 for Ai ^ 0. The solid line is a gaussian fit with a FWHM-width 
of 811 kHz. The peak Rabi frequencies are fii ^ 27r x 3 MHz and ^2 ~ 27r x 6 MHz. The 
error bars refer to the 1-sigma error in determining the particle number. 
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Fig. 4. A Ramsey-type experiment. Tlie population in the initial state \a> oscillates as 
the hold time Th during which both transfer lasers are off is increased. The solid line is a 
sinussoidal fit to the data up to r/i = 20 /is. Its frequency / is 115(2) kHz. The thin lines 
are borders to a region that is given by varying / by ±6 kHz, illustrating the estimated 
jitter in the two-photon detuning IA2 — Ai|. B Comparison of the rate of expansion in the 
horizontal direction for the molecular sample without and with STIRAP transfer. The top 
curve (circles) shows the Thomas-Fermi radius r of the molecular sample as a function of 
expansion time without STIRAP. The linear fit gives a rate of expansion of dr/dt = 1.0(1) 
mm/s, corresponding to an energy of ks x 14(4) nK. The bottom curve (squares) shows 
the expansion after double STIRAP with dr/dt = 0.7(1) mm/s, corresponding to ks x 7(2) 
nK. 
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